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Abstract 

The coefficients of different mass formulae derived from the liquid drop model and 
including or not the curvature energy, the diffuseness correction to the Coulomb 
energy, the charge exchange correction term, different forms of the Wigner term 
and different powers of the relative neutron excess / = (N — Z)/A have been de- 
termined by a least square fitting procedure to 2027 experimental atomic masses. 
The Coulomb diffuseness correction Z 2 j A term or the charge exchange correction 
Z 4 / 3 /A 1//3 term plays the main role to improve the accuracy of the mass formula. 
The Wigner term and the curvature energy can also be used separately for the same 
purpose. The introduction of an \1\ dependence in the surface and volume energies 
improves slightly the efficiency of the expansion and is more effective than an / 
dependence. Different expressions reproducing the experimental nuclear charge ra- 
dius are provided. The different fits lead to a surface energy coefficient of around 
17-18 MeV and a relative equivalent rms charge radius ro of 1.22-1.23 fm. 
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1 Introduction 



The binding energy of new nuclides both in the superheavy element 
region and the regions close to the proton and neutron drip lines are still 
poorly known and the different theoretical extrapolations do not fully agree. 
Therefore other studies are still necessary to better predict the masses of such 
exotic nuclei. The atomic nucleus resembling to a charged liquid drop, semi- 
macroscopic models including a pairing energy term have been firstly advanced 
to determine the nuclear masses [Tf2"] . Macroscopic- microscopic approaches, 
mainly the finite-range droplet model and the finite-range liquid drop model 
[3] have been proposed to simulate the smooth part of the nuclear masses 
and the non smooth part depending on the parity of the proton and neutron 
numbers and the proximity of the magic numbers. Nuclear masses have also 
been reproduced accurately within the statistical Thomas-Fermi model within 
a Seyler-Blanchard effective interaction [^|5] . Microscopic Hartree-Fock self- 
consistent calculations using the mean-field approach and Gogny or Skyrme 
forces and pairing correlations [6|7] as well as relativistic mean field theories 
[S] have also been advanced to try to reproduce these nuclear masses. Neural 
networks have been recently used [9J in the same purpose. 

Beyond the description of the nuclear ground state energy, the evolution 
of the nuclear binding energy with deformation and rotation governs the fis- 
sion, fusion, cluster and a decay potential barriers and the existence of the well 
deformed rotating states. Within the liquid drop model approach, the main 
characteristics of these barriers may be reproduced using, firstly, four basic 
macroscopic terms : the volume, surface, Coulomb and nuclear proximity en- 
ergy terms and, secondly, shell and pairing energy contributions to explain 
structure effects and improve quantitatively the results |1U|I11II12|I13|I14|I15] . 

In a first paper [16J some coefficients and terms of the liquid drop model 
mass formula and a specific nuclear radius have been fitted using a set of 1522 
experimental nuclear masses. The purpose of the present work is to extend 
this investigation to determine the relative efficiency of different more com- 
plex combinations of terms of the liquid drop model to reproduce the masses 
[TT] of 2027 nuclei and to study, particularly, the separated influence of the 
curvature energy, the different forms of the Wigner term and of different pow- 
ers of the relative neutron excess I = (N — Z)/A. As suggested in [18] fits on 
ground state masses alone have been prefered as was done with the Hartree- 
Fock mass formulae to avoid that the parameters might be distorted by the 
theoretical and experimental uncertainties associated with the barriers. The 
second aim is to fit the experimental nuclear charge radii by different expres- 
sions and to compare with the relative charge radii tq = Rq/A 1 ^ derived 
from the different mass formulae. Finally, the coefficients of the mass formu- 
lae using these different expressions for the charge radius are determined as 
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well as their accuracy. Another motivation of this work is to improve later 
the coefficients of the generalized liquid drop model previously proposed 
Recently, the mutual influence of terms in semi-empirical formulae has been 
deeply investigated [T9] . 



2 Nuclear binding energy 



The nuclear binding energy B nuc i(A,Z) is the energy needed for separat- 
ing all the nucleons constituting a nucleus. It is related to the nuclear mass 
M n , m by 

B nud (A, Z) = Zmp + Nm N - M n . m (A, Z). (1) 

B nMC /(A,Z) may thus be connected to the experimental atomic masses given 
in [IT] since : 

M n . m (A, Z) = M a . m (A, Z) - Zm e + B e (Z). (2) 
The binding energy B e (Z) of all removed electrons is [18] 



B e (Z) =a el Z 2 - 39 + b el Z 5 - 35 , (3) 

with a d = 1.44381 x 1(T 5 MeV and b el = 1.55468 x lCT 12 MeV. 

Different subsets of the following expansion of the nuclear binding energy 
in powers of A^ 1 ^ 3 and |/| have been considered : 

B = a v (l - k Vl \I\ - k V2 I 2 - k V3 I 4 ) A - a s (l - fc ai |/| - k s J 2 - k s J 4 ) A§ 
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3 e 2 Z 2 Z 2 Z% 
-a k (l - k kl \I\ - k k J 2 - hj 4 ") A^ - a A° - -— h f p — + a c , exc —Y 

V ' Kq A A 3 

Epair E s h e n E\\rig ner . (4) 

The first term gives the volume energy corresponding to the saturated ex- 
change force and infinite nuclear matter. I 2 A is the asymmetry energy of the 
Bethe-Weizsacker mass formula. The second term is the surface energy. It takes 
into account the deficit of binding energy of the nucleons at the nuclear surface 
and corresponds to semi-infinite nuclear matter. The Bethe-Weizsacker mass 
formula does not consider the dependence of the surface energy on /. This 
term was originally contained in the Weizsacker formula pQ. The third term is 
the curvature energy. It is a correction to the surface energy resulting from lo- 
cal properties and consequently depending on the mean local curvature. This 
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term is considered in the Lublin-Strasbourg Drop (LSD) model [20J, the TF 
model [5] but not in the FRLDM [3] . In the three first terms a dependence on 
|/| and I 4 , the so-called malacodermous term, has been envisaged since they 
have been proposed to better reproduce the fission barrier heights [2TJ and to 
simulate the softening of the surface of highly neutron-rich nuclei [22]. The 
A term appears when the surface term is extended to include higher order 
terms in A~ 1//3 and /. The fifth term gives the decrease of binding energy due 
to the repulsion between the protons. In the Bethe-Weizsacker mass formula 
this Coulomb energy is proportional to Z(Z — 1). Different formulae will be 
assumed for the charge radius. The Z 2 /A term is the diffuseness correction 
to the basic sharp radius Coulomb energy term (called also the proton form- 
factor correction to the Coulomb energy in [3]) and the term proportional to 
Z^'^/A 1 ' 3 is the charge exchange correction term. 

The pairing energy has been calculated with the following expressions 
used for spherical nuclei in the recent version of the Thomas- Fermi model [5]. 
For odd Z, odd N and N=Z nuclei 



E Pair = 4.8/N 1/3 + A.8/Z 1/3 - 6.6/ A 2/3 - 30/ A MeV. (5) 
For odd Z, odd N and N ^ Z nuclei 

E Pair = 4.8/N 1/3 + A.8/Z 1/3 - 6.6/ A 2/3 MeV. (6) 
For odd Z, even N nuclei 

E Pair = A.8/Z 1 ' 3 MeV. (7) 
For even Z, odd N nuclei 

E Pair = A.8/N 1 ' 3 MeV. (8) 
For even Z, even N nuclei 

Ep air = 0. (9) 

The theoretical shell effects used in the TF model (7 th column of the 
table in [4] and [3]) have also been retained since they allow to reproduce 
correctly the masses from fermium to Z = 112 [23]. They have been calcu- 
lated from the Strutinsky shell-correction method and previously to the other 
coefficients of the TF model. The fits on nuclear masses depend necessarily 
on the choice of the selected theoretical shell effects or the formulae chosen to 
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describe these shell effects. The sign for the shell energy term comes from the 
adopted definition in [4]. It gives, for example, a contribution of 12.84 MeV 
to the binding energy of 208 Pb. 

The Wigner energy allows to reproduce the kink in the nuclear mass 
surface that is not a shell effect in the usual sense. It depends on I and 
appears in the counting of identical pairs in a nucleus. Different expressions 
are considered. The first expression is simply W\I\ [24J. Its effect is to decrease 
the binding energy when N 7^ Z. 

The congruence energy term is given by : 



= -10exp(-4.2|/|) MeV. (10) 

It represents an extra binding energy associated with the presence of congruent 
pairs [5]. 



Within an Hartree-Fock approach [23] it has been assumed that there 
is nothing compelling about an exponential representation and a gaussian 
expression 

E = V W exp(-XI 2 ) (11) 
is just as acceptable. 



Another term has also been proposed in [25 



E = (3\N - Z\ exp[-(A/A ) 2 \. (12) 
We have also tested separately two other possible expressions 

E = /3\N — Z\ exp [—(A/Aq)] . (13) 

and 



E = ^TTW < 14 > 

The term of nuclear proximity energy does not appear in the binding energy 
of the ground state since it becomes effective only for necked shapes but not 
for slightly deformed ground states. 
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3 Coefficients of the mass formulae 



To obtain the coefficients of the different expansions by a least square 
fitting procedure, the masses of the 2027 nuclei verifying the two conditions 
: N and Z higher than 7 and the one standard deviation uncertainty on the 
mass lower than or equal to 150 keV [17] have been used. These restrictions 
are not employed by all investigators [T9|26] . The root-mean-square deviation 
a defined by 



^ = S [M Th - M Exp f 
n 

has been used to determine the relative efficiency of the different selected 
sets of terms since n ^> / where f is the number of fit parameters. A very 
efficient software has been used. The extraction of standard errors on the fit 
parameters seems very difficult but the errors are surely very weak. On the 
other hand, the values of the last decimals of a coefficient can be changed in 
counterbalancing by a change in the last decimals of another coefficient for 
almost the same rms deviation. 



In Table I, the improvement of the nuclear mass reproduction when ad- 
ditional contributions are added to the basic A, AI 2 , A 2 / 3 , A 2 ^ 3 I 2 , Z 2 /A 1 ^ 3 
terms is clearly displayed. The curvature energy is not taken into account. 
The introduction of the pairing term is obviously necessary. The introduction 
of a constant term improves slightly the adjustment and changes strongly the 
surface energy coefficient. It induces also a severe discontinuity during the 
transition from one to two-body shapes as in fission, fusion or a emission. 
The congruence energy term at least with the fixed coefficients adopted here 
(as in the LSD and TF models) is much less efficient to lower a than the 
Wigner term W|J|. When the coefficients before the exponential and the ex- 
ponent are free the congruence energy tends to the usual Wigner term since 
the coefficient before the exponential diminishes while the exponent increases. 
The diffuseness correction to the Coulomb energy and the standard Wigner 
term W|J| can be used separately to strongly lower a. The Wigner energy 
is approximately independent of the nuclear shape [21]. In a division into 2 
fragments, all with the same value of |/|, the form W|/| for the Wigner energy 
jumps at scission to 2 times its original value, which leads to a discontinuity 
of around 6.7 MeV of the potential energy at the scission point between the 
nascent fragments for 258 Fm as an example. For the congruence energy the 
discontinuity is also important : 3.9 MeV for this same nucleus. The Coulomb 
diffuseness correction term has the main advantage to be almost continuous 
at the scission point in the entrance or exit channels. The combination of the 
Coulomb diffuseness correction term and the W|/| term allows to reach the 
very satisfactory value a = 0.608 [3T71I2T)] . The yl 2 / 3 |/| term is useful to improve 
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the accuracy of the expansion and is more effective than the A 2 / 3 I 4 term but 
when the Wigner term and the Coulomb diffuseness correction factor term are 
taken into account the introduction of the A 2 / 3 |/|, A 2 / 3 1 4 and A terms are 
ineffective. 

In Table II, the efficiency of the curvature energy term with different / 
dependences is investigated without taking into account the Wigner contribu- 
tion. When the Coulomb diffuseness correction is disregarded the introduction 
of the term in A 1 / 3 allows to decrease a of 0.1 MeV. When the Coulomb dif- 
fuseness correction is considered the term in A 1 / 3 is ineffective . The addition 
of a A 1 / 3 / 2 term improves slightly the results. Supplementary terms in |/| 
to calculate the volume, surface and curvature energies allow finally to reach 
a=0.58 MeV. They are still more efficient than the J 4 terms. The curvature 
energy has the advantage to be continuous at the scission point at least in 
symmetric entrance or decay channels. It has the disadvantage that its value 
(and its sign) lacks of stability. 

These two first tables show a good convergency of the volume a^ and 
asymmetry volume k v constants respectively towards around 15.5 MeV and 
1.7—1.9 . The variation of the surface coefficient is larger but a s evolves around 
17-18 MeV. Small values of the surface coefficient favors quasi-molecular or 
two-body shapes at the saddle-point of the potential barriers while large values 
of a s promote one-body elongated shapes. As it is well known the surface 
asymmetry coefficient k s is less easy to precise. 

For the Bethe-Weizsacker formula the fitting procedure leads to 



with a=1.30 MeV. That gives r =1.216 fm and k„=1.506. The a value is 
explained by the non dependence of the surface energy term on the relative 
neutron excess /. 

The mass formulae including a Wigner term given by the formulae 
(12), (13) or (14) are examined in Table III. The values of A which minimize 
the mass rms deviation are respectively 48, 35 and 40 for these three expres- 
sions. They are determined with an accuracy of around 5 mass numbers. These 
formulas for the Wigner energy are supposed to be approximately independent 
of the nuclear shape. Their discontinuity at the scission point of fission or fu- 
sion barriers is less important than that of the congruence and W|/| terms. For 
example, for the expressions \N — Z\ exp [— (A/48) 2 ], 1.5\N — Z\ exp [— (A/35)} 
and 1.6 ^^^y and symmetric decay of the 254 Fm the discontinuities at the 
contact point of the nascent fragments are respectively 0, 0.05 and 2.2 MeV. 



-0.71056 




B nucl (A, Z) = 15.5704A - 17.1215A 2/3 
~ 23.4496/ 2 A - E pair - E shell 



(16) 
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The first, fifth and eighth lines of the table indicate that within the sim- 
plest form for the volume and surface energies , thus with only 7 parameters, 
the r.m.s deviation on the masses is less than 0.6 MeV. The introduction of 
a dependence on |/| and of the curvature energy lowers a of 0.04 MeV within 
the expression (12) and of 0.02 MeV for the two last ones. The surface energy 
coefficients are generally lower with these forms of the Wigner energy and, 
consequently, is generally higher. 

Within the Wigner term derived from the Hartree-Fock approach |6J and 
taking A = 28 and A = 485 in the formulas (11) and (12) as derived using 
the BSk2 Skyrme force the following formula is obtained 



B = 15.4503 (l - 1.7463/ 2 ) A - 17.5701 (l - 1.5296/ 2 ) A 2/3 



0.6e 2 Z 2 Z 2 48 -, 

+ 1.1948— - E pair - E shell + 1.769e~ 485/ 
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1.219AV3 1 A pair 

-0.2197|V-Z|e- (A/28)2 (17) 

with a = 0.623 MeV. 

The accuracy can be improved in varying A and A : 



B = 15.4121 (l - 1.7972/ 2 ) A - 17.3059 (l - 1.7911/ 2 ) A 2/3 

O.Qe 2 Z 2 Z 2 mT 2 

- 232A i /3 + 0.8961— - E pair - E shell + 2.25e- 80/ 



-0.4883|iV- Z\e- {A/50) 

It leads to a — 0.573 MeV and produces sizable changes in several of the 
other parameters which proves the importance of the Wigner term and the 
mutual influence of terms in the semi-empirical mass formulae |19j . 



In Table IV the efficiency of the charge exchange correction term in 
Z 4/3/ A i/3 is s t u died. The lines can be compared respectively to the 7 th line of 
table 1, 2 nd line of table 2, 8 th line of table 1, 3 rd line of table 2, Q th and 12 th 
lines of table 1 and l si and 5 th lines of table 3. The only change is that the 
diffuseness correction term in Z 2 j A has been replaced by the charge exchange 
correction term. It is quite interesting to observe that the introduction of 
these two terms separately leads to the same accuracy, the same value of the 
charge radius, almost the same values of the surface coefficients and to small 
changes of the volume coefficients. So the charge exchange correction term 
is as efficient as the diffuseness correction term and there is no correlation 
between these two terms and the charge radius and the surface coefficients 
and weak correlation with the volume coefficients. The introduction of both 
these two terms does not allow to improve the accuracy of the mass formulae 
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and leads to spurious values of the volume coefficient. 

Finally, the shell and pairing energies, the diffuseness correction term 
to the Coulomb energy, the charge exchange correction term and the Wigner 
energy separately lead to significant reductions of the rms deviation while 
the other terms increase the accuracy by only 2 or 3 per cent. A correlation 
between the surface coefficient and the radius can be extracted. The relative 
radius r diminishes when the surface coefficient increases. 



4 Nuclear charge radius 

Experimentally and for nuclei verifying N and Z higher than 7 the set 
of 782 ground state nuclear charge radii presented in ref. [27J indicates a rms 
charge radii of 0.94944A 1 / 3 which leads for the equivalent rms charge radius 
(denoted by Q in ref. [28]) given by 



to the value i?o — 1-2257 A 1 / 3 fm with a = 0.124 fm. Other data have been 
obtained recently [29]. In the adjustment to the nuclear masses displayed in the 
tables I and II the reduced charge radius r charge converges to 1.22-1.23 fm 
in good agreement with the experimental data for the charge radius. The 
introduction of the expressions (12), (13) and (14) for the Wigner term leads 
to slightly higher values. 

The experimental data indicate that the ratio Rq/A 1 / 3 is not strictly 
constant. For example, Rq/A 1 ' 3 = 1.312 fm for 40 Ca and Rq/A 1 ' 3 = 1.234 
for 48 Ca while Rq/A 1 ' 3 = 1.217 for 190 Pb and Rq/A 1 /' 3 = 1.201 fm for 214 Pb. 
These Ro and r values correspond to the Coulomb energy of a charged liq- 
uid drop with constant charge density and spherical sharp surface. It must 
not be confused with half-density radius derived from two-parameter Fermi 
function or other parameters entering in other parametrisations of the charge 
distribution. 

In the adjustment to the experimental nuclear masses the nuclear mass 
radius is not fitted. Root-mean-squared matter radii are given in Ref. J3U] 
for specific nuclei. For all isotopic series a decrease of the mass rms radius 
is observed with increasing neutron number as for the charge radius. In this 
mass range (A = 63 — 75) the radius of the neutron and proton distributions 
are very similar. Global fits lead to an overall contraction of the nuclear radius 
as T = |iV — Z\/2 increases [21], a so called "isospin shrinkage". Recently, an 




(19) 
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estimate of the neutron skin linear in the relative neutron excess I has been 
extracted from experimental proton radii and observed mirror displacement 
energies |32j and the coefficients have been also determined from antiprotonic 
X-rays and radiochemical data [35] . 

Within a simple form of the mass formula including the Coulomb dif- 
fuseness correction and the Wigner term given by the formula (13) but not 
the dependence on the curvature energy and on |/| and I 4 the introduction of 
the factor r = 1.2257 deduced from the experimental charge radius leads to 



B = 15.3543 (l - 1.7445/ 2 ) A - 17.2293 (l - 1.5765/ 2 ) A 2/3 

6e 2 Z 2 Z 2 

— + 1.2U2—-E pair -E shell -0.76U\N-Z\e- A/35 (20) 



1.2257AV3 ' A 
and a = 0.615 MeV. 

More accurate formulae may be used to reproduce the mean behaviour 
of the equivalent rms charge radius. 

The expression 



#0 = 1.0996A 1/3 + 0.653 fm (21) 

leads to o = 0.066 fm. Its introduction in the above-mentioned mass formula 
curiously diminishes the accuracy to a = 1.16 MeV. The same behaviour 
is observed when the dependence in Z 2 / A is replaced by a dependence in 
Z^/A l '\ 

The formula 



R = 1.1718A^ + ^^- fm (22) 

leads to o = 0.064 fm, while for the corresponding mass formula a = 0.929 MeV. 

The adjustment to the experimental charge radii within the form pro- 
posed in Ref. [34J gives 



R = 1.1818A 1/3 - 0.089 + L ^ 3 8 fm (23) 

with a = 0.064 fm and a = 0.900 MeV while the form selected in Ref. [35] 
leads to 
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^ = , 176M V3 + ^ + I^ /m (24) 

and a = 0.064 fm for the radius and 0.88 MeV for the corresponding mass 
formula. 

The form of the expression given the equivalent rms charge radius pro- 
posed in Ref. [36] 



#0 = 1.2332A 1 / 3 + —|— - 0.18688A 1 / 3 / fm (25) 
gives a good accuracy o = 0.052 fm but o = 0.765 MeV. 

The fact that improved formulae for the charge radius lead to signifi- 
cantly poorer fits to the masses calls into question the reliability of the equiv- 
alent uniform charge distribution as a link between charge radii and Coulomb 
energy. 

Finally a rms deviation for the charge radius of only 0.01 fm (fitted on 
362 nuclei) has been obtained in Ref. [31] in expressing the charge mean square 
radius by 



r 3 = a + bA + cT (26) 

for spherical nuclei and taking into account the deformation for the other 
nuclei via a nucleonic promiscuity form factor P depending on the distance of 
Z and N to the proton and neutron magic numbers. 

To describe the main properties of the fusion, fission, cluster and a emis- 
sion potential barriers in the quasi-molecular shape path, the formula 



#0 = 1.28A 1/3 - 0.76 + 0.8A- 1/S (27) 

proposed in Ref. [M] for the equivalent rms radius has been retained in the 
GLDM [H)lfll1fT2l^^ This formula coming from the Droplet Model 

approach simulated a small decrease of the density with increasing mass. This 
does not seem corroborated by the experimental data [27] at least on the charge 
radius. Nevertheless without changing the surface and Coulomb energies which 
are shape- dependent terms, the adjustment of the volume energy coefficients 
and the introduction of a pure Wigner term leads to 



B = 15.5209 ( 1 - 1.93AI 2 ) A - 17.9439 ( 1 - 2.QI 2 ) A 2/3 



0.6e 2 Z 2 Z 2 



1.28AV3 _ 0.76 + 0.8A-V3 ■ A 



+ 2.1798— - E pair - E shell - 27.21|/| (28) 
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and a = 0.686 MeV. The introduction of other terms can still allow to dimin- 
ish a. 



The Fig. 1 displays the difference between the theoretical and experi- 
mental masses within a formula of Table I recalled below : 



B 



0.6e 2 Z 2 

O28lW + L0339 A 



15.3543 (l + 0.02841/1 - 1.8837/ 2 ) A 
17.0068 (l + 0.1311/1 - 2.274Z 2 ) A 2/3 
Z 2 



— E pair — E she u — 16.27|/| 



(29) 



with o = 0.605 MeV. For most of the nuclei with A higher than 110 the 
difference between the theoretical and experimental masses is less than 1 MeV. 




Fig. 1. Difference (in MeV) between the theoretical and experimental masses for the 
2027 selected nuclei as a function of the mass number. 



5 Summary and conclusion 



The coefficients of different mass formulae derived from the liquid drop 
model and including or not the curvature energy, the Coulomb diffuseness 
correction energy, the charge exchange correction term, different forms of the 
Wigner term and different powers of the relative neutron excess / have been 
determined by a least square fitting procedure to 2027 experimental atomic 
masses. 

The Coulomb diffuseness correction Z 2 /A term or the charge exchange 
correction Z 4 ^ 3 /A 1 ^ 3 term plays the main role to improve the accuracy of the 
mass formula. The Wigner term and, with a smaller effect, the curvature en- 
ergy can also be used separately for the same purpose but their coefficients are 
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unstable. The introduction of an \I\ dependence in the surface and volume en- 
ergies improves only slightly the efficiency of the expansion. An J 4 dependence 
and a constant term are not worth including. Expressions for the Wigner term 
leading to small discontinuities at the scission point of the potential barriers 
are as efficient as the W|J| term. Different expressions reproducing the exper- 
imental nuclear charge radius are provided. The different fits lead to a surface 
energy coefficient of around 17-18 MeV and a relative equivalent rms charge 
radius r of 1.22-1.23 fm. 

References 

[I] C.F. von Weizsacker, Z. Physik 96 (1935) 431. 

[2] H.A. Bethe, R.F. Bacher, Rev. Mod. Phys. 8 (1936) 82. 

[3] P. Moller, J.R. Nix, W.D. Myers, W.J. Swiatecki, At. Data Nucl. Data Tables 
59 (1995) 185. 

[4] W.D. Myers, W.J. Swiatecki, LBL report 36803, 1994. 

[5] W.D. Myers, W.J. Swiatecki, Nucl. Phys. A 601 (1996) 141. 

[6] M. Samyn, S. Goriely, P.-H. Heenen, J.M. Pearson, F. Tondeur, Nucl. Phys. A 
700 (2002) 142. 

[7] J. Rikovska Stone, J. Phys. G 31 (2005) R211. 

[8] M. Bender, et al., Phys. Lett. B 515 (2001) 42. 

[9] S. Athanassopoulos, E. Mavrommatis, K.A. Gernoth, J.W. Clark, Nucl. Phys. 
A 743 (2004) 222. 

[10] G. Royer, B. Remaud, J. Phys. G 10 (1984) 1057. 

[II] G. Royer, B. Remaud, Nucl. Phys. A 444 (1985) 477. 
[12] G. Royer, J. Phys. G 26 (2000) 1149. 

[13] G. Royer, R. Moustabchir, Nucl. Phys. A 683 (2001) 182. 

[14] R.A. Gherghescu, G. Royer, Phys. Rev. C 68 (2003) 014315. 

[15] C. Bonilla, G. Royer, Heavy Ion Phys. 25 (2006) 11. 

[16] G. Royer, C. Gautier, Phys. Rev. C 73 (2006) 067302. 

[17] G. Audi, A.H. Wapstra, C. Thibault, Nucl. Phys. A 729 (2003) 337. 

[18] D. Lunney, J.M. Pearson, C. Thibault, Rev. Mod. Phys. 75 (2003) 1021. 

[19] M. W. Kirson, Nucl. Phys. A 798 (2008) 29. 



13 



[20] K. Pomorski, J. Dudek, Phys. Rev. C 67 (2003) 044316. 

[21] M. Dahlinger, D. Vermeulen, K.H. Schmidt, Nucl. Phys. A 376 (1982) 94. 

[22] A.K. Dutta, J. P. Arcoragi, J.M. Pearson, R.H. Behrman, M. Farine, Nucl. Phys. 
A 454 (1986) 374. 

[23] S. Hofmann, et al., Z. Phys. A 354 (1996) 229. 

[24] W.D. Myers, Droplet Model of Atomic Nuclei, Plenum, New York, 1977. 

[25] S. Goriely, M. Samyn, P.-H. Heenen, J.M. Pearson, F. Tondeur, Phys. Rev. C 
66 (2002) 024326. 

[26] A.E.L. Dieperink, P. Van Isacker , Eur. Phys. J. A 32 (2007) 11. 

[27] I. Angeli, At. Data Nucl. Data Tables 87 (2004) 185. 

[28] R.W. Hasse, W.D. Myers, Geometrical relationships of macroscopic nuclear 
physics, Springer- Verlag, Berlin, 1988. 

[29] J. Libert, B. Roussiere, J. Sauvage, Nucl. Phys. A 786 (2007) 47. 

[30] G.F. Lima, et al., Nucl. Phys. A 735 (2004) 303. 

[31] J. Duflo, Nucl. Phys. A 576 (1994) 29. 

[32] J. Duflo, A.P Zucker, Phys. Rev. C 66 (2002) 051304(R). 

[33] E. Friedman, A. Gal, J. Mares, Nucl. Phys. A 761 (2005) 283. 

[34] J. Blocki, J. Randrup, W.J. Swiatecki, C.F. Tsang, Ann. of. Phys. 105 (1977) 
427. 

[35] P.A. Seeger, W.M. Howard, Nucl. Phys. A 238 (1975) 491. 
[36] W.D. Myers, W.J. Swiatecki, Phys. Rev. C 62 (2000) 044610. 



14 



Table 1 

Dependence of the energy coefficient values (in MeV or fm) on the selected term set 
including or not the pairing and congruence terms and root mean square deviation. 
The theoretical shell energies are taken into account but not the curvature energy. 
The Coulomb energy is determined by a c -rrm with a c = 3e 2 /5rQ. 



a v 


k Vl 


ky 2 


ky 3 


a s 


ks 1 


k S2 


k S3 a 


W 


Cong 


Pair 


^0 


fp 


a 


15.5959 


- 


1.7051 




17.1723 


- 


0.994 






n 


n 


1.227 




1.32 


15.4996 


- 


1.7011 




16.7628 


- 


0.990 






n 


y 


1.235 




0.977 


15.0775 


- 


1.6799 




14.8548 


- 


1.005 


7.686 




n 


y 


1.269 




0.917 


15.93 


- 


1.8454 




18.7399 


- 


1.722 


- 




y 


y 


1.204 




0.870 


15.3416 


- 


1.8220 




16.0795 


- 


1.838 


10.982 




y 


y 


1.251 




0.728 


15.8559 


- 


1.8549 




19.3863 




1.999 






y 


y 


1.191 


1.1478 


0.628 


15.4108 




1.7119 




17.5361 




1.404 






n 


y 


1.218 


1.3755 


0.660 


15.4721 




1.7154 




17.8591 




1.411 


-1.2 




n 


y 


1.212 


1.4358 


0.659 


15.4436 




1.5909 


2.147 


17.4801 




0.732 


10.96 




n 


y 


1.215 


1.2286 


0.640 


15.4367 


-0.0551 


1.9205 




17.2772 


-0.296 


2.420 






n 


y 


1.215 


1.1095 


0.616 


15.3306 




1.8991 




16.1189 




2.410 




41.82 


n 


y 


1.254 




0.729 


15.351 




1.8097 




16.9773 




1.972 




21.60 


n 


y 


1.233 


0.9621 


0.608 


15.3543 


-0.0284 


1.8837 




17.0068 


-0.131 


2.274 




16.27 


n 


y 


1.228 


1.0339 


0.605 


15.5494 




1.8406 




17.9723 




2.077 


-4.130 


26.08 


n 


y 


1.216 


1.0681 


0.590 


15.3619 




1.7539 


0.799 


17.0142 




1.682 


3.88 


19.57 


n 


y 


1.230 


0.9879 


0.606 



Table 2 

Dependence of the energy coefficient values (in MeV or fm) on the selected term 
set and root mean square deviation. The theoretical pairing and shell energies are 
included but not the Wigner energy. The Coulomb energy is determined by 
with a c = 3e 2 /5ro. 



a v 


k Vl 


k V2 


k V3 a s k si 


k 82 


k S3 


afc 


k kl 


k k2 


kk 3 r 


fp 


a 


14.772 




1.6753 


12.3434 


1.094 




7.602 






1.285 




0.873 


15.4623 




1.7140 


17.8861 


1.402 




-0.564 






1.214 


1.4157 


0.660 


15.444 




1.8730 


18.0077 


2.783 




-1.269 




47.98 


1.221 


1.3721 


0.639 


15.4401 




2.3713 


-7.447 17.9965 


7.136 


-65.24 


-1.983 




134.19 


-1530.7 1.230 


1.1196 


0.594 


15.7016 


0.1622 


1.351 


19.6025 1.223 


-1.445 




-5.116 


11.917 


-27.66 


1.220 


1.0166 


0.583 



Table 3 

Dependence of the energy coefficient values (in MeV or fm) on the selected term 
set and root mean square deviation. The theoretical pairing and shell energies are 
included. The Wigner energy determined by the expressions (12), (13) or (14) are 
taken into account. The I A term is disregarded. 







k V 2 




ks\ 




"'ft 


ku 


ku 


(12) 


(13) 


(14) 


' u 


fn 
JP 


(j 


15.2565 




1.7475 


16.8485 




1.603 








0.679 






1.241 


1.1161 


0.600 


15.311 


-0.035 


1.8745 


16.826 


-0.197 


2.238 








0.535 






1.235 


0.9642 


0.580 


15.5431 




1.6662 


18.8035 




0.721 


-3.343 




-12.95 


1.017 






1.220 


1.2826 


0.580 


15.8011 


0.0752 


1.5133 


20.0578 


0.466 


-0.223 


-5.754 


4.068 


-15.75 


0.794 






1.213 


0.9836 


0.556 


15.2374 




1.7880 


16.6394 




1.803 










1.514 




1.245 


0.9639 


0.590 


15.2949 


-0.0047 


1.8077 


16.7004 


-0.057 


1.958 










1.319 




1.241 


0.8456 


0.585 


15.3474 


0.1082 


1.3103 


17.3787 


0.800 


-2.111 


-1.789 


16.015 


-85.97 




2.249 




1.247 


0.7978 


0.571 


15.2738 




1.8238 


16.6658 




2.010 












1.262 


1.242 


0.8822 


0.581 


15.315 




1.7455 


16.7943 




1.150 


-0.217 




-190.59 






1.581 


1.238 


0.8382 


0.574 


15.265 


0.1212 


1.2866 


16.8666 


0.873 


-2.401 


-1.041 


22.288 


-154.51 






2.516 


1.254 


0.7380 


0.563 



Table 4 

Dependence of the energy coefficient values (in MeV or fm) on the selected term 
set and root mean square deviation. The theoretical pairing and shell energies are 
included as well as the charge exchange correction term. The Coulomb energy is 
determined by a c -^^ with a c = 3e 2 /5ro. 



a v 


k V2 


a s 


k S2 




kk 2 
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W 


Cong 


0(12) 


0(13) 


ro 


&c,exc 


a 


15.2393 


1.7168 


17.5359 


1.404 
















1.218 


1.2987 


0.660 


15.2876 


1.7191 


17.8984 


1.402 


-0.584 














1.214 


1.3361 


0.660 


15.2949 


1.7206 


17.8647 


1.411 






-1.222 










1.212 


1.3530 


0.659 


15.2752 


1.8799 


18.0208 


2.785 


-1.291 


47.24 












1.221 


1.2934 


0.638 


15.7124 


1.8602 


19.3863 


1.999- 










y 






1.191 


1.0850 


0.628 


15.2305 


1.8138 


16.9776 


1.972 








21.58 








1.233 


0.9100 


0.608 


15.1186 


1.7518 


16.8486 


1.603 












0.678 




1.241 


1.0501 


0.600 


15.1174 


1.7910 


16.6395 


1.803 














1.514 


1.245 


0.9082 


0.589 



